High-pT Tomography oi d + Au and Au + Au at SPS, RHIC, and LHC 
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The interplay of nuclear effects on the pr > 2 GeV inclusive hadron spectra in d + Au and 
Au + Au reactions at \/sj^j^ = 17, 200, 5500 GeV is compared to leading order perturbative QCD 
calculations for elementary p + p{p + p) collisions. The competition between nuclear shadowing, 
Cronin effect, and jet energy loss due to medium-induced gluon radiation is predicted to lead to a 
striking energy dependence of the nuclear suppression/enhancement pattern in A + A reactions. We 
show that future d + Au data can used to disentangle the initial and final state effects. 

PACS numbers: 12.38. Mh; 24.85. +p; 25.75.-q 



Introduction. Tomography is the study of the proper- 
ties of matter through the attenuation pattern of fast par- 
ticles that propagate and lose energy as a result of multi- 
ple elastic and inelastic scatterings. Recently, this tech- 
nique has been applied in the field of nuclear physics [1,2] 
to map out the evolution of the QCD matter density pro- 
duced in ultra-relativistic heavy ion reaction. It is based 
on the theoretical advances in understanding QCD multi- 
parton interactions in non-Abelian media [3,4]. 

The determination of the opacity of the transient 
quark-gluon plasma (QGP) produced in such reactions 
via jet tomography requires theoretical control over the 
interplay between many competing nuclear effects that 
modify the high-pT hadron spectra. These include nu- 
clear modifications to the parton distribution functions 
(PDFs), referred to as shadowing [5], Cronin effect [6], 
and jet quenching [7], as well as the energy dependence 
of the underlying pQCD parton spectra. In this letter we 
propose an approach to disentangle these effects by com- 
paring high-pT hadron yields in p-l-p (p-|-p), d+A (p+A), 
and A + A reactions over a very wide energy range. In 
particular, predictions are presented for d + Au and cen- 
tral Au + Au at center of mass energies per nucleon 
Vsattv = 17, 200, 5500 GeV typical of the CERN Super 
Proton Synchrotron (SPS), the Relativistic Heavy Ion 
Collider (RHIC), and the future Large Hadron Collider 
(LHC). 

We demonstrate that at SPS the puzzling absence of 
quenching of tt" in central Pb + Pb [8] can be under- 
stood as due to a larger than previously expected Cronin 
enhancement [9] dominating over our predicted [4] sup- 
pression due to jet energy loss. At RHIC energies, on 
the other hand, we find that quenching dominates over 
both Cronin and shadowing effects. Furthermore, the 
interplay of these effects leads to a surprising approxi- 
mately constant suppression factor of the Glauber geom- 
etry scaled [10] pQCD prediction in the 4 < pT < 20 GeV 
range. At LHC we predict that the 7r° suppression fac- 
tor is substantially larger than at RHIC but also de- 
creases systematically with transverse momentum in the 
6 < Pt < 100 GeV range due to the higher initial gluon 
densities expected and the hardening of the underlying 
initial jet spectra. 



Particle spectra in d + A and A + A. The scaling of 
high-pr hadron production in d + A and A + A is simply 
controlled by nuclear geometry in the absence of initial 
and final state interactions. The Glauber multiple col- 
lision model [10] can be used to calculate the number 
of binary nucleon-nucleon collisions at any impact pa- 
rameter b. In p + A the experimental cross section has 
usually been presented without centrality selection while 
in A -|- A reactions hadron multiplicity distributions are 
generally presented with restricted centrality (impact pa- 
rameter b) cuts. Dynamical nuclear effects for these cases 
are detectable through the nuclear modification ratio 
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where 2 A and Taa(&) = / d'^r Ta{y)Tb{y - b) in terms 
of nuclear thickness functions Ta(^) = J dz pyi(r, z) are 
the corresponding Glauber scaling factors of daPP. The 
lowest order pQCD differential cross section for inclusive 
A + B ^ h + X production that enters Eq.(l) is given 
by 
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In Eq.(2) Xa,Xb are the initial momentum fractions car- 
ried by the hard-scattered partons with probabilities 
sampled from the PDFs fa/A{xa,Qa)- The momen- 
tum fraction carried away by the leading hadron Zc = 
Ph/Pc is sampled from the fragmentation functions (FFs) 
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Dh/c{zc,Qc)- We use the leading order (LO) Gliick- 
Reya-Vogt (GRV98) parameterization of PDFs [11] and 
the LO Binnewies-Kniehl-Kramer (BKK) parameteriza- 
tion of the FFs [12]. 

A constant if-factor and parton kx broadening func- 
tion, g(k) — exp(— kp/(k|n)pp)/7r(k|n)pp, are included 
to account phenomenologically for next-to-leading order 
corrections. In systematic fits [13] to the inclusive hadron 
production iiip+p{p+p) reactions these parameters are 
fixed from data [6,14]. The ii'-factor drops out in the 
ratios oi B + A to scaled p + p, but the phcnomenolog- 
ical kx broadening in p -\- p is essential to establish an 
accurate nuclear geometry scaled baseline. We find that 
a fixed {k'^)pp = 1.8 GeV^ reproduces to within 30% the 
spectral shapes in p + p{p + p) for px > 2 GeV over the 
whole energy range of interest. 

In B + A reactions isospin effects are accounted for on 
average in the PDFs for a nucleus with Z protons and 
iV neutrons via fa/Aixa,Ql) = {Z/A) fa/p{xa,Ql) + 
(N/A) fa/n{xa,Qa)- Nuclcar modifications to the 
PDFs [5] are included through the shadowing function 
SA{xa,Qa) from thc EKS98 parameterization [15]. 

The Cronin effect observed in p -I- ^ reactions relative 
to the Glauber-scaled p + p result [6] is modeled via mul- 
tiple initial state scatterings of the partons in cold nuclei. 
For an initial state parton distribution dN^^\\i), random 
elastic scattering induces further /c^-broadening as shown 
for example in [9]. The possibility of hard fluctuations 
along the projectile path leads to a power law tail of the 
kx distribution that enhances (Afc|,)-^ beyond the naive 
Gaussian random walk result x^^, where x = {n) — L/X 
is the cold nuclear opacity in terms of the path length 
L and the parton mean free path A. The screening scale 
/i regulates the Rutherford divergence in a cold nucleus. 
For XM^ < Q'^ax the Rutherford tail leads to a 

logarithmic enhancement of the mean square momentum 
transfer {/S.k?^)^ — xAi^lii(l -I- cQ^^^/^^), where c de- 
pends on the detailed form of the kinematic cut-off. For 
a high energy parton with transverse momentum px pro- 
duced in a, p + A reaction Q^^^ ^ Px- We therefore 
model the Cronin effect by using 
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{kx)pA ~ {f^x)pp 
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in the kx broadening functions g{ka) in Eq.(2), tak- 
ing La = 4:/3Ra as the mean nuclear thickness tra- 
versed. Fig. 1 shows that the calculation is consistent 
with the energy dependence y/s = 27 A, 38.8 GeV [6] 
observed in p + W/p + Be with transport parameters 
set as follows: c/^u^ = 0.18 /GeV^ A = 3.5 fm and 
/\ = 0.05 GeV^/fm. These are consistent with ^j? /\ — 
0.064 ± 0.036 GeV^/fm extracted from fits to Drell-Yan 
data [16]. Fig. 1 also shows that the expected Cronin -|- 
shadowing effect in p -I- W/p -f Be at RHIC energies is 
much smaller than at lower energies because the high-p^ 
pQCD spectra at y/sjqpf = 200 GeV are considerably less 
steep. 
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FIG. 1. The ratio of y4-scaled p + W/p + Be data on tt^, tv~ 
production at ^/s — 27.4,38.8 GeV from [6]. Calculations 
for ^(tt^ + tt^) include nuclear shadowing and initial parton 
broadening as in Eq.(3) with ^i^/A = 0.05 GeV^/fm. The 
anticipated y/s — 200 GeV p + W/p + Be ratio is also shown. 



The effects of multiple scattering and nuclear shadow- 
ing in d -|- Au and Au + Au without final state interac- 
tions at SPS, RHIC, and LHC are shown in Fig. 2 for 
neutral pions. The numerical results for charged parti- 
cles are comparable. Variations arise from the different 
partonic contribution and the correspondingly different 
shadowing for various hadron species [13]. In our model 
of thc Cronin effect, the enhancement in Au + Au at SPS 
energies of \/s^^ = 17 GeV may reach a factor ~ 4 
aX px ~ 4 — 5 GeV. This is greater than observed in 
Ph + Ph reactions and also greater than estimated with 
the Cronin model of Ref. [8]. Unfortunately, at these 
low energies the results are extremely sensitive to model 
assumptions due the very rapid fall-off of the partonic 
spectra. We note that at least within our model, there 
is room for hadron suppression due to energy loss even 
at SPS. At RHIC for ^Sj^j^ = 200 GeV the Cronin en- 
hancement spans the pt = 1 ^ 8 GeV range and is seen 
to peak at px — 3 GeV. Its maximum value in d + Au 
and Au -t- Au is 1.3 (1.6) respectively. Similar and even 
smaller magnitudes of the Cronin effect at RHIC have 
been discussed in [17]. At higher transverse momenta 
the effects of isospin and shadowing lead to Rba — 0.8 
at Px ^ 20 GeV. At LHC energies of Vsatjv = 5500 GeV 
Cronin effect is overwhelmed by shadowing at small x 
(px < 10 GeV) and anti-shadowing at larger x when the 
EKS98 [15] parameterization is used. The net nuclear 
modification due to Cronin effect and shadowing at LHC 
is expected to be tiny (< 15%) throughout the px range 
shown. 

We turn next to the predicted suppression effects in 
nucleus-nucleus reactions due to jet quenching [7]. In 
Eq.(2) this is taken into account in the fragmentation 
function via the modification of thc momentum fraction 
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carried away by the leading hadron. If a jet of momen- 
tum pc prior to hadronization looses a fraction < e < 1 
of its energy then z — pu/pc z* — z/(l — e). The 
distribution P(e, E) of the fractional energy loss of a fast 
parton with energy E due to multiple gluon emission is 
computed as in Gyulassy-Levai-Vitev (GLV) [2]. 
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FIG. 2. The nuclear modification Rba{pt) due to Cronin 
efi'ect and shadowing (but not energy loss) for tt" in d + Au 
{B = d,A = Au) and central Au + Au {B = A = Au) reac- 
tions at Vsjvjv = 17, 200, 5500 GeV. 



We compute P(e, E) taking into account the longi- 
tudinal Bjorken expansion the plasma (gluon) density 
p(t) = (to/t)p(to), where tqPo {l/nR\)dNS /dy re- 
lates to the gluon rapidity density produced in cen- 
tral A + A that fixes the initial opacity. It has been 
shown that the azimuthally averaged energy loss is in- 
sensitive to transverse expansion [18]. The mean num- 
ber of radiated gluons {N3{E)) remains small due to 
the plasmon mass cut-off uipi ^ 0.5 GeV [2]. There- 
fore, there is a finite n = (no radiation) contribution 
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6{e). We have checked the sensitiv- 



ity of the results to reducing the plasmon mass by a factor 
of two. This was found to lead to ^ 25% more suppres- 
sion at Pt = 5 GeV and to < 10% increased suppression 
at PT = 20 GeV for RHIC energies. 

Our main results for central Au+Au including all three 
nuclear effects (Cronin-f Shadowing-|-Quenching) are pre- 
sented in Fig. 3. Jet tomography consists of determining 
the effective initial gluon rapidity density dN^/dy that 
best reproduces the quenching pattern of the data [8,19- 
21]. At SPS the large Cronin enhancement is reduced 
by a factor of two with dN^ /dy = 350 but the data are 
more consistent with a smaller gluon density < 200. Un- 
fortunately, as emphasized above, at this low energy the 
results are very sensitive to the details of the model. At 
RHIC, ioi Pt > 2 GeV jet quenching dominates, but sur- 
prisingly the rate of variation with pT of the Cronin en- 
hancement and jet quenching conspire to yield an approx- 
imately constant suppression pattern with magnitude de- 
pendent only on the initial dN^ /dy. At higher pT > 20 
GeV the softening of the initial jet spectra due to the 
EMC modification of the PDFs compensates for the re- 
duced energy loss. This unexpected interplay between 
the three nuclear effects at RHIC is the main prediction 
of this letter. At LHC energies the much larger gluon 
densities dN^/dy ^ 2000 — 3500 are expected to lead to 
a dramatic variation of quenching with pT as shown. 

In nuclear media of high opacity the mean fractional 
energy loss ( Ai?) / E oi moderately hard partons can be- 
come on the order of unity. For LHC this may be reflected 
in the pt < 10 GeV region through deviations from the 
extrapolated high-p^ suppression trend. Hadronic frag- 
ments coming from energetic jet would tend to compen- 
sate the rapidly increasing quenching (seen in Fig. 3) with 
decreasing transverse momentum and may restore the 
hydrodynamic-like participant scaling in the soft regime. 

Conclusions. In this letter we predicted a character- 
istic evolution pattern of the magnitude and the pT de- 
pendence of the nuclear modification factor in d + A and 
A + A reactions as a function of the center of mass en- 
ergy per nucleon. A systematic approach was used to 
take into account Cronin effect, nuclear shadowing, as 
well as jet quenching. Our results suggest that at SPS 
energies the Cronin enhancement may be larger than ex- 
pected previously, leaving room for moderate energy loss. 
At RHIC we predict that the three nuclear effects in cen- 
tral Au -\- Au lead to a surprising approximately constant 
suppression pattern of tt" with RAAipr) — 0.3 — 0.2 
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for dNs/dy ~ 800 - 1200. We emphasize that none 
of the nuclear effects alone would lead to such a flat 
Raa{pt)- At LHC shadowing and Cronin effect in the 
6 < Pt < 100 GeV range were found to be essentially 
negligible, leading to < 15% correction, while the jet 
quenching was predicted to be large and with a strong 
Pt dependence. We emphasize the importance of future 
d + Au data at RHIC to isolate and test the initial state 
Cronin and shadowing effects predicted in Fig. 2. While 
it is still too early to draw conclusions from the prelim- 
inary data [20,21] shown in Fig. 3, the combined future 
analysis of c? + Au and Au + Au high-pT measurements 
will improve the tomographic determination of the initial 
gluon densities produced at RHIC. 
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FIG. 3. The suppression/enhancement ratio Raa{pt) 
{A — B — Au) for neutral pions at \/sjyjv ~ ^"^^ 200, 
5500 GeV. Solid (dashed) lines correspond to the smaller 
(larger) effective initial gluon rapidity densities at given ^/s 
that drive parton energy loss. Data on tt" production in cen- 
tral Pb + Pb at VSiviv = 17.4 GeV from WA98 [8] and in 
central Au + Au at v^jyjv ~ ^^O GeV [19], as well as pre- 
liminary data at 200 GeV [20] from PHENIX and h"^ cen- 
tral/peripheral data from STAR [21] are shown. The sum of 
estimated statistical and systematic errors are indicated. 
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